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Abstract 

If entanglement is available, the error-correcting ability of quantum codes can be 
increased. We show how to optimize the minimum distance of an entanglement-assisted 
quantum error-correcting (EAQEC) code, obtained by adding ebits to a standard 
quantum error-correcting code, over different encoding operators. By this encoding 
optimization procedure, we found several new EAQEC codes, including a family of 
[[n, 1, n; n — 1]] EAQEC codes for n odd and code parameters [[7, 1, 5; 2]], [[7, 1, 5; 3]], 
[[9, 1, 7; 4]], [[9, 1, 7; 5]], which saturate the quantum singleton bound for EAQEC codes. 
A random search algorithm for the encoding optimization procedure is also proposed. 
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1 Introduction 



Since Shor proposed the first quantum error-correcting code [T], the theory of quantum error 
correction has been extensively developed. Today, quantum stabilizer codes [21 El IH El E] are 
the most widely-used class of quantum error-correcting codes. One reason for this is that 
the CSS and CRSS code constructions [71 El El E] allow classical self-orthogonal codes to be 
easily transformed into quantum stabilizer codes. 

Bowen constructed the first entanglement-assisted quantum error-correcting (EAQEC) 
code from a three-qubit bit-flip code with the help of two pairs of maximally-entangled states 
[U] . Bowen's code, which can correct an arbitrary one-qubit error, serves as an example that 
entanglement increases the error-correcting ability of quantum error-correcting codes. Brun, 
Devetak and Hsieh showed that if shared entanglement between the encoder and decoder is 
available, classical linear quaternary (and binary) codes that are not self-orthogonal can be 
transformed to EAQEC codes [101 E]. 

An [n, k, d] classical linear quaternary code encodes k quaternary information digits into 
n quaternary digits and can correct up to [^"J quaternary digit errors, where d is called the 
minimum distance of the code. Brun, Devetak and Hsieh showed that an [n, k, d] classical 
linear quaternary code can be transformed to an [[n, 2k — n + c, d; c]] EAQEC code that 
encodes Ik — n + c information qubits into n qubits with the help of c pairs of maximally- 
entangled states (ebits) for some c [TO]. This EAQEC code can correct up to |_^r"J qubit 
errors and has the same minimum distance d as the classical code. If entanglement is used, 
it boosts the rate of the code. However, it has not been explored how entanglement can 
instead help increase the minimum distance. In addition, given parameters n, k, c, it is not 
clear how to construct an [[n, k, d; c]] EAQEC code directly. We will answer these questions 
in the paper. 

An EAQEC code can be obtained from a standard QECC by changing one or more an- 
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cilia qubits of the initial state with ebits. We first discuss how adding the maximum number 
of ebits introduces symplectic partners of the generators of the stabilizer group, and define a 
selection operator that determines the set of logical operators. Then the minimum distance 
of the EAQEC codes can be optimized over distinct selection operators. If we add fewer than 
the maximum ebits, we have the freedom to choose the set of generators of the stabilizer 
group, and the freedom to replace different ancilla qubits with ebits. EAQEC codes can be 
optimized over these choices. These factors, together with the selection operator, can be 
represented in terms of unitary row operations. Applying this encoding optimization proce- 
dure to some standard quantum stabilizer codes, we construct several new EAQEC codes, 
including a family of [[n, l,n;n — 1]] EAQEC codes for n odd, and codes with parameters 
[[7, 1, 5; 2]], [[7, 1, 5; 3]], [[9, 1, 7; 4]], [[9, 1, 7; 5]], which saturate the quantum singleton bound 
for EAQEC codes, and are not equivalent to any standard quantum stabilizer code. (We say 
that an [[n, k, d; c]] EAQEC code is not equivalent to any standard quantum stabilizer code 
if there is no standard [[n + c, k, d]] quantum code.) We also compare the EAQEC codes 
obtained by the encoding optimization procedure with the EAQEC codes obtained by the 
construction of [TO] . 

Although the encoding optimization procedure seems to be a promising method to 
construct EAQEC codes with high minimum distance d for given parameters n, k, c, it will 
be shown that the complexity of the encoding optimization procedure increases exponentially 
with n + k, which implies that it is impossible to fully optimize the minimum distance for 
high n + k. Hence we develop a random search algorithm, which achieves a suboptimal 
result efficiently On the other hand, we perform a different computer search for EAQEC 
codes that have a circulant check matrix and find several EAQEC codes that achieve the 
quantum singleton bound and are not equivalent to any standard quantum stabilizer code. 
This circulant construction of EAQEC codes also provides evidence that entanglement helps 
increase the error- correcting ability of the quantum error-correcting codes. 

This paper is organized as follows. Basics of stabilizer codes and EAQEC codes are 
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introduced in Section [2j We discuss the encoding optimization procedure by first considering 
the case of maximal entanglement and the selection operator in Section 13.11 including the 
construction of the family of [[n, l,n;n— 1]] EAQEC codes for n odd. We then generalize 
to adding arbitrary amounts of entanglement. The effect of unitary row operations are 
discussed in Section 13. 2^ which completes the encoding optimization procedure. The results 
of applying the encoding optimization procedure to some standard quantum stabilizer codes 
are provided in Section 13.31 Then we propose the random search algorithm for EAQEC 
codes in Section HJ together with several examples. The circulant construction of EAQEC 
codes is described in Section EJ followed by the discussion section. 

2 Preliminaries 

2.1 Stabilizer Codes 

Suppose S is an abelian subgroup of Pauli group Q n that does not include —I, with a set 
of r = n — k independent generators {gi, <?2, • ■ ■ , g r }- An [[n, k, d]] quantum stabilizer code 
C(S) corresponding to the stabilizer group S is the 2 k - dimensional subspace of the n— qubit 
state space fixed by S. The minimum distance d is the minimum weight of an element in 
J\f(S) — S, where A/ r («S) is the normalizer group of S. An element g = i m M\ ® M2 <8> • • • <8> M n 
in Q n , where Mj G {/, X, Y, Z} and m £ {0, 1, 2, 3}, can be expressed as g = i m> X a Zp with 
a,/3 two binary n-tuples and m' £ {0,1,2,3}. In this expression, if Mj = I, X, Z, or Y, 
then the j—th. bits of a and (3 are (atj,f3j) = (0,0), (1,0), (0, 1), or (1, 1), respectively, and 
m! = m + l (mod 4), where I is the number of Mj's equal to Y. We define a homomorphism 

¥ >(? m X a Zp) = (a,f3), 
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and define a symplectic inner product between two elements and (a-}, fit) in Z; 

by 



2/i 



(«!, ft) © (a 2 , /? 2 ) fo + Pi - "2, 

where • is the usual inner product in Z£ . Thus, two elements g, h in (? n commute if and only 
if the symplectic inner product <£>(g) © p(h>) is zero. Then a check matrix if corresponding 
to the stabilizer S is defined as a binary r x 2n matrix such that the the i— th row vector of 
if is (fi{gi)- For convenience, if is denoted by 

where Hx and Hz are two r x n binary matrices. The check matrix H must satisfy the 
following commutative condition, 

HA 2n H T = H x Hz + H z Hx = O rxr , 



where A 



2d 



Ifixn O n xn 

identity matrix. The i-th co 



O 



nxn ± nxn 



Oi X j is an % x j zero matrix, and i rxr is an r-dimensional 



umn of Hx is the error syndrome of the error operator Zj and 
the j-th column of Hz is the error syndrome of the error operator Xj. The error syndrome 
of Y\ is the sum of the /-th column of Hx and the /-th column of Hz- In general, the 
error syndrome of g = i m Mi ® M 2 <8> ■ ■ ■ <E> M n in Q n , where Mj G {i, X, Y, Z}, is the linear 
combination of the error syndromes of each factor Mj. For a code with minimum distance 
d, if the error syndromes of error operators of weight smaller than or equal to |_^"J are 
distinct, we call that code nondegenerate. Otherwise, it is degenerate. 
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The encoding procedure is described as follows. Consider the initial n— qubit state 



|0> |0> - - - 10> 



r=n— k 



where there are r = n — k ancilla qubits |0)'s and an arbitrary A;— qubit state 
generators of the stabilizer group of this class of states is 



>. A set of 



Z x =ZIIII ■ ■ ■ I 
Z 2 =IZIII 



Z r =I--IZI_^J / 

k 



with a check matrix 



The logical operators on are 



Or 



O r x(n—r) 



(1) 



and 



If Ue is a unitary operator such that {UeZiU e , • ■ ■ , UEZ r U E } is a set of generators of the 
stabilizer group S, then Ue is an encoding operation of C(S), and the encoded state Ue 
is fixed by the stabilizer group S. In particular, we can choose 



gi = U E ZiW E 
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for % = 1, • • • , r. The logical operators on Ue are 



Zi — UEZ r+1 U E , ■ ■ ■ , Zk — UEZ n U E , 

X\ = UEX r+ iU E , ■ • • , Xk = \JeX t JJ e . 

Note that the logical operators commute with the stabilizers, and the normalizer group of S 
is 

N [S] = (<?i, #2, • ■ ■ > Or-, Zi, Z 2 , • ■ • , z k , X l5 x 2 , • ■ ■ , 

with dimension In — r = r + 2fc. 

Given a check matrix if of a stabilizer group, Wilde gave an algorithm [12] to find an 
encoding circuit for this quantum stabilizer code. This algorithm applies a series of CNOT 
gates, Hadamard gates, Phase gates, SWAP gates, and row operations to the check matrix 
H such that H takes the form ([T]). This process is like performing Gaussian elimination on a 
matrix, but using CNOT gates, Hadamard gates, Phase gates, and SWAP gates, in addition 
to the elementary row operations of Gaussian elimination. There are two types of elementary 
row operations over the binary field: adding one row to another, which corresponds to mul- 
tiplying an operator by another; and exchanging two rows, which corresponds to relabeling 
two generators. Performing these row operations does not change the row space and hence 
the error-correcting ability of the codes. The effects of these gate operations on the entries 
in the check matrix follows: 

1. A CNOT gate from qubit i to qubit j adds column i to column j in Hx and adds 
column j to column % in Hz. 

2. A Hadamard gate on qubit % swaps column % in Hz with column % in Hx- 

3. A Phase gate on qubit i adds column i in Hx to column i in Hz- 

4. Three CNOT gates implement a SWAP gate. The effect of a SWAP gate on qubits % 
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and j is to swap columns i and j in both Hx and Hz- use the check matrices such 
that no Row operations in the encoding algorithm. 

The series of operations used in the algorithm serve as a unitary operation U E such that 
UeQiUe — an d hence the inverse operator \J E is a desired encoding operation. 

Note that Ue is not unique. The encoding algorithm generates, in general, different 
U E 's for different sets of generators of the same stabilizer group; and performing the steps 
in the encoding algorithm in a different order will, in general, generate a different U e even 
for the same set of generators. 

2.2 Entanglement- Assisted Quantum Error-Correcting Codes 

Brun, Devetak and Hsieh proposed a theory of quantum stabilizer codes when shared entan- 
glement between the encoder (Alice) and decoder (Bob) is available in [TO] . Suppose that 
Alice and Bob share c pairs of maximally entangled states. Assume further that Bob's halves 
of the c ebits are not subject to any error. Now if an arbitrary subset T = {t 1; • ■ ■ , t c } of the 
r = n — k ancilla qubits in are replaced by these c pairs of maximally-entangled states 
|$ + } j4B, s, then the (n + c)-qubit initial state is 



For convenience, the qubits on Alice's side will be numbered 1 to n and the qubits on Bob's 
side will be numbered 1 to c. Hence the tj-th qubit of Alice and the i-th qubit of Bob form 
a pair of maximally-entangled state. Then a set of independent generators of a stabilizer 



r 



\1>)ea = ® fa) ® \<P) , 



where 




S 



group of \ip) EA is 

I B , if i i T; 

for i = 1, • • • , r, 

<g> Zf , if i = tj G T; (2) 

X*®X?, forj = l,---,c. 

Note that the operators on the left and right of the tensor product are applied to Alice's 
qubits and Bob's qubits, respectively, and the superscripts A and B will be omitted through- 
out the rest of this article. The logical operators on \ip) EA are Z r+ i <S> I, • ■ ■ , Z n <g> I, and 
X r+1 J, • • • , X n <S> I. Now consider the operators on Alice's qubits. It can be observed 
that 

[Zi,Zj\ = ti, forO<i,j<r, (3) 
[X ti ,X tj ] = 0, forO<i,j<c, (4) 
{Z ti ,X u } = 0, forO <i<c, (5) 
[Z*,^.] = 0, forz^j, (6) 

where [g, h]= gh — hg and {g, h} = gh + hg, or 

^(Zi) </?(^i) = 0, for < i, j < r, (7) 

<p{X u ) <^(X f .) = 0, for < i, j < c, (8) 

tp(Z ti )0<p(X ti ) = l, for0<*<c, (9) 

<p(Zi) Q <p(X t .) = 0, for i^V (10) 

If a set of (r + c) operators satisfy equations ©-© or equations (17|)- (|TU|) . we say that the 
two operators in (JSj) or the two vectors in (Q form a symplectic pair, and they are symplectic 
partners of each other. Hence Z t . and form a symplectic pair. 

An encoding operation Ue is applied on Alice's n qubits, while no operation is performed 
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on Bob's c qubits. A set of generators of a stabilizer group S of the encoded state (Ue ® 
J ) \^)ea is i9ii ■ ■ ■ ,gr,h, ■ ■ ■ A}, where 

U E ZiUl ® I, if z i T; 
U E ZiU E ®Zj, if i = tj G T 

and 

h j = \J E X t] U E ®X j . 
The logical operators on {Ue® I) I^ea are 

Z x = U E Z r+1 U E ® /,-■-, Z k = U E Z n U E ® /, 

and 

X, = U E X r+1 U E ® /,..., X fc = f/ B X re 4 ® /. 

The 2 fc — dimensional subspace of the (n + c)— qubit state space fixed by the stabilizer 
group S with a set of generators {g±, ■ ■ ■ ,g r ,hi,--- , h c } is called an entanglement-assisted 
quantum error-correcting (EAQEC) code with parameters [[n, k, d; c]] for some minimum 
distance d. With the help of c pairs of maximally-entangled states, the stabilizer group 
of an [[n, k, d; c]] EAQEC code has c more generators than that of an [[n, k, d]] standard 
quantum stabilizer code. The operators on Bob's side have no direct effect on the error- 
correcting ability of the EAQEC code if we assume that the c qubits of Bob suffer no error. 
These operators serve to make the stabilizer abelian, so that the generators can be measured 
simultaneously. For convenience, we denote 

g[ = U E ZiUl 
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and 

h! 3 = XJ E X t U E , 

and the g[s and h'jS will be called the simplified generators of the stabilizer group. Similarly, 
we denote Z\ = UEZ r+i U E ,X'j = UEX r+ jU E . It is obvious that {g[, ■ ■ ■ ,g r ,h[,--- ,h' c } 
satisfy the commutative relations ((31)-®, and g' t . and are a symplectic pair. Let S' = 
{g'11 • " " ) g' r i ■ ■ ■ , h' c ), and Sj = (gj : j T). Hence the normalizer group of S' is 

M(S') = (g l : l ^T,Z[ r -- ,44- ,X' k ) 

with dimension 2n — (r + c) = 2/c + r — c. The minimum distance d of the EAQEC code 
defined by S is the minimum weight of an element in A/"(«S') — Sj. In particular, when c = r, 
we have Sj = and 

Af(S') = (Z[,--- ,Z' k ,X[,--- ,X' k ). 
An [[n, fc, d; c]] EAQEC code must satisfy the quantum singleton bound for EAQEC codes 

pro]: 

n + c — k> 2(d - 1). (11) 



We define a simplified check matrix H' as a binary (r + c) x 2n matrix such that the 
r + c row vectors of H' are ip(g'i) for i = 1, • • • , r and <^(/ij) for j = 1, • ■ • , c. For simplicity, 
we usually order the generators g[ and h'j so that <p{g'i) is the i— th row vector of if' for 
i = 1, ■ • • , r, (f(hj) is the (j + r)— th row vector of if' for j = 1, • ■ ■ , c, and the j-th and 
(j + r)— th row vectors are a symplectic pair. H' must satisfy the commutation relations 
(j71)-(fTUl). and in the case c = r, 



H'K 2n H' T 



O r x r 


irxr 


irxr 


O r xr 



(12) 
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For example, the simplified check matrix corresponding to the set of generators (j2J) of a 
stabilizer group of the initial state \ip) EA is 



O r xn 


I r xr Orx(n—r) 


I rxr O rx (n—r) 


Orxn 



(13) 



Conversely, an (r + c) x 2n binary matrix H, serving as a simplified check matrix, can define 
a stabilizer group and hence an EAQEC code. The number of ebits required to construct an 
EAQEC code p3] is 

c=-iank(HAH T ). (14) 
2 

Like a check matrix of a standard quantum error-correcting code, the simplified check matrix 
H' can be used to determine the minimum distance of nondegenerate EAQEC codes. Note 
that Wilde's encoding circuit algorithm [12] can also be applied to a simplified check matrix 
to find an encoding unitary operator of the EAQEC code, just as for a standard stabilizer 
code. 

Similarly, we define a simplified logical matrix L' corresponding to the logical operators 
by putting f{Z^) to be the %— th row vector of V for i = 1, • • • , k, and <p(X'j) to be the 
(j + k)—th row vector of V for j = 1, • • • , k. Since the logical operators commute with 

Wi, ■ ■ ■ ,g' r ,Ki ■ ■ ■ X)> we have 

H'A 2n L' T = O. (15) 
Since the logical operators satisfy the commutative relations ([3])-©, we have 



Okxk 


Ikxk 


Ikxk 


Okxk 
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For example, the simplified logical matrix corresponding to the initial state \ip) E ^ is 



Okxn 


Okxr Ikxk 


Okxr Ikxk 


Okxn 



(16) 



3 The Encoding Optimization Procedure for EAQECCs 

A standard [[n, 2k — n, d}} quantum stabilizer code, obtained by the CSS construction [TJ IE] 
from an [n, k, d] classical linear self-orthogonal code, has the same minimum distance as 
the classical code used. On the other hand, an [[n, 2k + c — n, d; c}} EAQEC code can be 
constructed from an [n, k, d] classical linear quaternary code by the construction of |10j . and 
c is determined by (TH|) . It seems that only the number of information qubits is increased 
by introducing ebits. However, with the help of entanglement it is possible to define more 
distinct error syndromes for a given codeword size, and hence the set of correctable error 
operators might be larger. We would like to construct EAQEC codes with a higher minimum 
distance instead of higher rate. 

One way to construct an EAQEC code is to start with a standard QECC and move c 
of the qubits from Alice's side to Bob's side. So long as c < d/2, the resulting code can be 
encoded by a unitary operator on Alice's side, given c ebits of initial shared entanglement 
between Alice and Bob. While such codes can be interesting and useful, they are not the 
subject of interest for this paper; because such codes retain an ability to correct errors on 
Bob's qubits, they are in a sense not making full use of the fact that Bob's halves of the 
ebits are noise-free. They therefore may not have the maximum error correcting power on 
Alice's qubits for the given parameters n, k and c. We are interested in EAQEC codes that 
can do better than any standard code in this sense. 

To make this idea precise, we say that an [[n, k, d; c]] EAQEC code is not equivalent to 
any standard quantum stabilizer code if there is no standard [[n + c, k, d}} quantum code. 
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If there exists a standard [[n + c, k, d}} quantum code, then we may not be achieving the 
maximum boost to our error correcting power from the c ebits of shared entanglement. We 
expect added entanglement in general to increase the error-correcting ability of a quantum 
error-correcting code such that the EAQEC code is not equivalent to any standard quantum 
stabilizer code, and indeed this turns out to be possible by our encoding optimization pro- 
cedure. (Note that this is not always possible — the smallest examples of the [[3, 1, 3; 2]] and 
[[4, 1, 3; 1]] codes are both equivalent to the standard [[5,1,3]] QECC, and this is the best 
that can be done.) 

We now consider how added entanglement affects standard quantum stabilizer codes. 
Suppose {g'^g'2, • • • ,g' r } forms a set of independent generators of a stabilizer group S for an 
[[n, k,d}} quantum stabilizer code C(S), and Ue is the encoding unitary operator obtained 
by Wilde's encoding circuit algorithm [12], such that UeZ\U e = g[,-" iUeZ t U e = g' r . 
The logical operators are Z[ = U E Z r+1 U E , ■■■ , Z' k = U E Z n U E , X[ = U E X r+1 U E , ■■■ , 
X' k = UeX h U e . Suppose entanglement between the encoder and decoder is available, and a 
set T of c ancilla qubits are replaced by ebits. An EAQEC code is obtained. This introduces 
c simplified generators h'j = U E X t .U E 's, tj 6 T, j — 1, • ■ ■ , c, to the generating set of the 
stabilizer group S. As we will examine in detail below, the encoding unitary operator is 
not uniquely defined. The EAQEC code defined by S' = (g[,--- , , g' r , h[, • ■ ■ , , h' c ) may gain 
higher error-correcting ability by modifying the encoding operator. 

We first discuss the case c = r, where the generators h'^s are symplectic partners of 
g-'s, respectively. We will treat the case c < r later, by optimizing the choice of c linearly 
independent generators from the group {h[, ■ ■ ■ , , h' r ). 

3.1 Selecting Symplectic Partners and Logical Operators 

Since the symplectic partners of g[, ■ ■ ■ ,g' r are not unique, we now explain how to select 
these partners such that the minimum distance of the EAQEC codes is higher than the 
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code without entanglement. Suppose W is a unitary Clifford operator that commutes with 
Zi, ■ ■ ■ , Z r such that after the operation of W, the simplified check matrix of the initial state 
( TT3l) becomes 



O r xn 


I r xr O r x{n—r) 


I rxr A 


C B 



(17) 



where A and B are two r x (n — r) binary matrices, and C is an r x r binary matrix. The 
simplified check matrix satisfies the commutation relations (1T])- (ll0j) with 

C T + AB T + C + BA T = O. (18) 

In addition, it can be checked that the simplified logical matrix is of the form 



Okxn 


A T Ikxk 


Okxn Ikxk 


B T Okxk 



after Gaussian elimination such that ( II 5 p and (TTUj) hold. Since 

(U E W)Z i (U E WY = U E ZiUl, = g[ 

for i = 1, • • • ,r, UeW is also an encoding operator of the quantum stabilizer code C(S). 
However, the symplectic partners of g^s, 

U E {WXiW ] )U E 

may differ from UEXiU E for % = 1, • • • , r and the logical operators 

U E (WXiWi)U E , U E {WZ^)U ] E , 
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for i,j = r+l,--',n are different. Choosing a set of matrices A, B, C such that C +AB + 
C + BA T = O determines a unitary operator W by the encoding circuit algorithm, which 
determines a set of symplectic partners of g[, ■ ■ ■ ,g' r and a different set of logical operators. 
Thus we call W the selection operator for EAQEC codes. The minimum distance of the 
EAQEC code can be optimized by examining each distinct encoding operator UeW. Note 
that the simplified logical matrix is not affected by the matrix C. Therefore, there are 2 2rk 
distinct sets of logical operators. 

Lemma 1. Given matrices A and B, then the matrix C, satisfying f llSp . is of the form 

C = BA T + M, 

where M is a symmetric matrix. 

Proof. Suppose C = BA T + M + M', where M 1 is not a symmetric matrix, satisfies (118j) . 
Then 

O = AB T + BA T + C + {C'Y = M' + (M') T , 

which implies that M' is the zero matrix. □ 

We construct an EAQEC code that achieves the quantum singleton bound by applying 
this procedure to a standard quantum stabilizer code in the following example. 

Example 1. A check matrix of the standard [[5, 1, 1]] 5-qubit bit flip code is 



00000 


11000 


00000 


01100 


00000 


00110 


00000 


00011 
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Applying the encoding circuit algorithm to this check matrix, we obtain an encoding operator 
Ue- In particular, if C = O in f iTHj) . then 



AB T + BA T = O. 



When k = 1, AB T + BA T = O holds if and only if A = B or at least one of A and B 
is the zero vector. Let W be the selection operator determined by the encoding circuit 



algorithm with A 







and B 
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Then the encoding operator 



UeW generates a [[5, 1, 5; 4]] EAQEC code with a simplified check matrix 



and a simplified logical matrix 



00000 


11000 


00000 


01100 


00000 


00110 


00000 


00011 


01111 


00000 


11000 


00000 


00011 


00000 


11110 


00000 


11111 


00000 


00000 


11111 



With the help of 4 ebits, the minimum distance is increased from 1 to 5. The quantum 
singleton bound ( fill) is saturated by the parameters [[5, 1,5; 4]]. Because the minimum 
distance of a standard [[9,1]] quantum stabilizer code is at most 3, this [[5, 1,5; 4]] code is 
not equivalent to any standard 9-qubit code. □ 



In Example [TJ we constructed a [[5, 1, 5; 4]] EAQEC code from a [[5, 1, 1]] 5-qubit bit flip 
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code, which is a quantum version of the [5, 1, 5] classical repetition code. The result can be 
generalized to the construction of an [[n, l,n;n — 1]] EAQEC code for n odd from an [n, l,n] 
classical repetition code as follows. 

Theorem 2. There are [[n, l,n;n— 1]] EAQEC codes for n odd that achieve the quantum 
singleton bound (TTTT) and are not equivalent to any standard quantum stabilizer code. 

Proof. Suppose H n is an (n — 1) x n parity-check matrix of a classical [n, l,n] repetition 
code: 

1 1 

1 1 '•• ••■ 

1 1 

1 1 

Then the [[n, 1, 1]] n— qubit bit-flip code has a check matrix 













We want to introduce (n — 1) simplified generators to the generating set of the stabilizer 
group such that the minimum distance of the code is increased to n. Consider a simplified 
check matrix 

O H n 

H' = 

H n O 

By (I14jl . the number of symplectic pairs in H' is 

^iank(H'AH' T ) = mnk(H n H^) = n - 1, 
for n odd. It can be verified that H' is a simplified check matrix with minimum distance n. 



H n = 
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Therefore, there exists a set of symplectic partners of the generators of the stabilizer group 
of the n— qubit bit flip code such that the minimum distance of the code is n. It is easy 
to verify that ( TTTT) is saturated by the parameters [[n,l,n;n — 1]]. These [[n,l,n;n — 1]] 
codes are not equivalent to any standard quantum stabilizer code, for there are no standard 
[[2n — 1, 1, n}} quantum codes. □ 

Example 2. By Theorem [2] with n = 5 and 



110 
110 
110 
1 1 



which is a parity-check matrix of a classical [5, 1,5] repetition code, a set of generators of 
the stabilizer of the [[5, 1, 1]] 5-qubit bit flip code is 

g[ =ZZIII 
g' 2 =IZZII 
g' 3 =IIZZI 
g' 4 =IIIZZ. 

Then we add the following set of generators: 

ti x =XXIII 
h 2 =IXXII 
h 3 =IIXXI 
ti 4 =IIIXX. 



These generators do not form symplectic pairs. Observe that [gi,gj] = 0, [hi,hj] = 0, and 

[g~i,hj] = except {gi,h} = 0, {g 2 ,hi} = 0, {g 2 ,h} = 0, {# 3 ,M = 0, {g 3 ,h 4 } = 0, 
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{g A ,h 3 } = 0. Let h[ = h 2 h A = IXXXX, h' 2 = h u h' 3 = h A , and h' 4 = hxh 3 = XXXXI. 
Then we have 



[g'i,9j] = 0, for all i, j, 



[//,/'] = 0, foralU,j, 



{^,//} = 0, foraUz, 
[g' t J'] = 0, forz^j. 



Hence we have introduced 4 simplified generators such that there are 4 symplectic pairs. 
Observe that if any one of the simplified generators h[, h ; 2 , h 3 , h' A is removed and c = 3, the 
minimum distance instantly drops to 2. If two simplified generators h[, h' 2 or h[, h' A or h' 2 , h' 3 



are removed and c = 2, the minimum distance further decreases to 1. 



□ 



According to [10], given a parity-check matrix H of an [n, k, d] classical binary linear 
code, an [[n, 2k + c — n,d; c]] EAQEC code can be constructed from a simplified check matrix 
H' defined as 



H' 



O H 
H O 



(19) 



where the number of ebits c required for this EAQEC code is given by ffT4l) . The family of 
EAQEC codes in Theorem [2] can also be obtained by this construction. When c = n — k, 
the quantum singleton bound (TTTT) becomes 



n — k > d — 1, 



which is exactly the same as the classical singleton bound. However, there are no nontrivial 
classical binary codes achieve the singleton bound from 
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3.2 Unitary Row Operators 



Since we have the freedom to choose among different sets of generators of a stabilizer group, 
and also the freedom to choose which ancilla qubits are replaced by ebits when c < r, we 
will show that the minimum distance can be further optimized over these two factors when 
c < r. We first discuss the effect of some "unitary row operators," which preserve the overall 
commutation relations (J3j)- (J6]) . 

Consider a unitary operator U = (I + iQ), where Q is a Pauli operator with eigen- 
values ±1. It is easy to verify that 



UgW 



9, if[Q,g] = 0; 
tQg, ii{Q,g} = 0. 



We define V\ 2 = V 3 V 2 Vi, where 



and 



Then 



t 9W2, ifj = 2; 
9-, ifjV2. 

Therefore, Vi^ is a unitary operator that performs multiplication of g[ to g' 2 , which corre- 
sponds to adding the first row to the second in the simplified check matrix. On the other 
hand, 

KK, ifj = l; 
h'j, if j ^1. 

Hence a row operation performed on {g[,--- ,g' r } induces a row operation performed on 
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Vi, 2 tijVl 2 



{hi, ■ ■ ■ , h' r } in order to preserve the commutation relations (JHI)- (J6j) . We call VJ.,2 a unitary 
row operator. Later we will need unitary row operators that perform the multiplication of gf i 
to h'j, the multiplication of Z[ to h'j, and the multiplication of X[ to h'j, respectively. These 
four types of unitary row operators are summarized in Table [TJ 



Table 1: Four types of unitary row operators 



Type 1. 


Vh'jV^ = 


/ h[h m , — h \r n i\r\ J 

I K, a j + i. V9jV ~\ 


g' m 9i, tij = m; 
g'j, if j + m. 


Type 2. 


Vh'jVi = 


I V~, l i vh'v^ = \ 

[ fy, if j + 1. i \ 


h' m 9i, if 3 = m; 
h'j, iij + m. 


Type 3. 


vhy = { 


h>, if j + l VX ^~- 


\ 9[X' m , if 3 = m- 
[X'j, iij + m. 


Type 4. 


Vh'^ = | 


\X , ifj = /; x 
K, if 3 + 1 VZ > V ~ 


J g'A, if j = m; 
1 Z'j, if j + m. 



When a different set of generators of the stabilizer group is chosen instead of {g[, ■ ■ ■ ,g' r }, 
this is equivalent to performing a unitary row operator V, which comprises a series of unitary 
row operators of type 1 on {g[, ■ ■ ■ ,g' r }. The operation of V on the simplified check matrix 
H' corresponding to {g[, ■ ■ ■ , g' r , h[, ■ ■ ■ , h' r } is to multiply H' by a (2n — 2k) x (2n — 2k) 
matrix of the form 



M z 








M x 



If 

M x = R m Rm-i •■•Ri, 
where i?£s are elementary row operations, then 

M z = R^Rm-i ' ' 'Ri- 

It can be checked that MH' satisfies ( fT2l) . If a set T = {t±, • • • , t c } of c < r ancilla qubits 
are replaced by ebits, it is possible that after the operation of V, the group S' x = (gj : j ^ T) 
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changes, and so does the set Af(S') — S' x . In addition, the span of a subset of {h[, • • • , h' r } 
can change after the operation of V, though the span of the full set remains unchanged. 
This means that if we add less than the maximum amount of entanglement to a code, we 
must optimize over such unitary row operations. Since the group S' x and the set N(S') — S' x 



for j ^ T, it suffices to assume that the operation V consists of type 1 unitary row operators 
that operate only on the h'j for j £ T. 

Let My be a c x r matrix such that the i— th row of My is the t— th row of Mz for 
i = 1, • • • , c. It is obvious that some My's have the same effect on the row space of H' . 
For example, if c = 2, {g[g' 2 , g' 2l - ■ ■ , g' r , h[, b! x b! 2 } and {g[, g' 2 , ■■■ ,g' r , h[, h' 2 } are two different 
sets of generators but they generate the same space and hence their corresponding EAQEC 
codes have the same minimum distance. Therefore, a distinct unitary row operation V is 
assumed to be be represented by a matrix My in reduced row echelon form. 

Theorem 3. The operation of V is equivalent to applying a series of type 1 unitary row 
operators on the hj for j £ T. In addition, there are 



l c =0 Z c _i=0 l c -2 h=0 

distinct unitary row operations. 

Proof. The total number of distinct unitary row operations N(r, c) is determined as follows. 
If we begin with matrices of the form 



remain the same under the operation of type 1 unitary row operators that operate on the hj 




1 



□ 



□ 



1 



□ 



□ 







1 □ 



□ 
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where □ can be or 1, there are 2 C ^'~ C ^ distinct unitary row operations. Now we consider 
matrices in which the leading ones are shifted to the right. Let lj denote the shift amount 
of the leading 1 of j-th row from its initial position for j — 1, • • • , c. It can be observed that 
lj < h if j < i- F° r a set {h, h,-'' > h}, the number of □ is c(r — c) — Yli=i lit ariG ^ hence 
there are 2 c ( r ~ c ^~^ i = 1 li distinct unitary row operations. Therefore, summing over all possible 
sets of {z'i, • • • , i c } shows that there is a total of 

r— c l c lc— 1 I2 

N(r,c) = Y: E E---E 2C(r " c) - ELl ' 

Z c =0Z c _i=0 l c -2 h=o 

distinct unitary row operations up to Gaussian elimination. □ 
Example 3. For r = 4 and c = 2, possible matrices of My are 



1 





□ 


□ 




1 


□ 





□ 




1 □ 


□ 










1 





□ 







1 


□ 













1 








1 


□ 


□ 










1 


□ 










1 










1 


□ 













1 













1 



There are 

2 4 + 2 3 + 2 2 + 2 2 + 2 1 + 2° = 35 
distinct row operations. □ 

The function N(r, c) has a symmetric property as in the following lemma and some 
closed forms of N(r, c) are listed in Table [21 

Lemma 4. N(r, c) = N(r, r — c) for any r and < c < r. 

Proof. We prove this lemma by mathematical induction. Assume N(l, c) = N(l, I — c) holds 
for / = 1, • ■ • , r and c = 0, • • • , /. It suffices to assume c<r — c or 2c <r. Consider the 
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following c x (r + 1 — c) matrix: 



1 ••• □ 
1 ••• □ 







1 □ 



c+1 



□ ••• □ 

□ ••• □ 



□ 



□ 



Dividing the matrix into two parts, we find that 



N(r + 1, c) = N(c + 1, c)N(r, c). 



Now consider the following matrix: 



1 
1 



□ 
□ 



• □ 
■ □ 



Similarly, we have 



■•• 1 □ ••• □ 







1 □ 



r+1— c 



c 



□ 



N(r + 1, r + 1 - c) = N(c + 1, l)N(r, r - c). 



From the above two equations, we have N(r + 1, c) = N(r + 1, r + 1 — c). By induction, we 
obtain the lemma. □ 



On the other hand, it can be observed that the selection operator W in the previous 
subsection can be decomposed as a series of unitary row operators of type 2, type 3, and 
type 4, for the matrix A determines a series of type 4 unitary row operators, the matrix B 
determines a series of type 3 unitary row operators, and the symmetric matrix M, satisfying 
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Table 2: Closed forms of N(r, c) for different c 



c 


N(r, c) 


or r 


1 


1 or r — 1 


2 r - 1 


2 or r — 2 


2 at— 1 nr- 1 i ! 
3^ Z 3 


3 or r — 3 


1 or-1 2^-2 , lor- 2 1 
21 ° 3^ 3 Z 21 


4 or r — 4 


4 1«r-2 1 or-2 , 'I at -2, 1 or-3 i 1 
315 U 21° " r 9 21 315 


5 or r — 5 


1 qor— 2 4 i «r— 3 _i_ 1 cr— 3 2 Ar—i _i_ 1 or— 4 1 
9765° 315 U 63° 63 315 9765 


6 or r — 6 


8 R/|r— 3 1 qor— 3 i 4 -r cr— 4 1 or— 4 i 1 ai — 5 1 or— 5 i 
615195 U 9765°" 1 945 U 441° 1 945 9765" 1 615195 



C = BA + M, determines a series of type 2 unitary row operators. Actually, unitary row 
operators of type 2 do not affect the set Af(S') — S' x or the error-correcting ability, and the 
symmetric matrix M can be eliminated. If a set T = {t±, ■ • ■ ,t c } of c < r ancilla qubits 
are replaced by ebits, it can be verified that Af(S') = (gj : j T, Z 1: • • ■ , Z k , Z\, ■ ■ ■ , Z k ) 
remains unchanged under the operation of type 3 and type 4 unitary row operators on the h'j 
for j ^ T. It suffices to assume that the operation W consists of type 3 and type 4 unitary 
row operators that operate only on the hi, for j G T. To sum up, we have the following 
theorem. 

Theorem 5. The operation of W is equivalent to applying a series of type 4 unitary row 
operators, followed by a series of type 3 unitary row operators, on the hj for j G T. In 
addition, there are 2 2ck distinct selection operators with 

C = BA T . 

Combining the effects of the unitary row operation V with the selection operator W in the 
previous section, we can optimize an encoding operation of the form U = VUeW over 

2 2cA W(r, c) 

possibilities. We call this the encoding optimization procedure for EAQEC codes. Note that 
we can find another unitary row operator W corresponding to W such that W'Ue and UeW 
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are equivalent encoding operators. While W operates on the raw stabilizer generators and 
logcial operators, W operates on the encoded stabilizer generators and logical operators. 
Hence, we can also solve the optimization problem for an operator of the form U = VW'U e 
(which is what we actually do in practice, combining VW into a single optimization). 

3.3 Results of the Encoding Optimization Procedure 

We apply the encoding optimization procedure to a [[7,1,3]] quantum BCH code [15j US] 
and Shor's [[9, 1,3]] code pQ and the results are shown in Table |3] and Table SJ where d opt is 
the minimum distance of the optimized EAQEC codes, d st d is the highest minimum distance 
of an [[n + c, k}} standard stabilizer code, and N opt is the number of encoding operators that 
give an EAQEC code with minimum distance d opt . 

Example 4. The check matrix of a standard [[7, 1,3]] quantum BCH code adopted in the 
encoding optimization procedure is 

0000000 1001011 
0000000 0101110 
0000000 0010111 
1001011 0000000 
1100101 0000000 
1011100 0000000 

As shown in Table El the parameters [[7, 1,7; 6]], [[7, 1,5; 3]] and [[7, 1,5; 2]] achieve the 
quantum singleton bound for EAQECC ( ITTT) and are not equivalent to any standard quantum 
stabilizer code. We would like to compare these two EAQEC codes to a competing EAQEC 
code with n = 7 and d = 5 by the construction of [10] or the binary version (fT9|) . According to 
Grassl's table p2], a classical linear code over GF(A) (or GF(2)) that meets our requirement 
is a [7, 2, 5] linear quaternary code, which can be used to construct a [[7, 2, 5; 5]] EAQEC 
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Table 3: Optimization over the [[7, 1, 3]] quantum BCH code 



c 


dopt 


dstd 


N op t 


2' 2ck N(r, c) 


6 


7 


5 


36 


4096 


5 


5 


5 


31920 


64512 


4 


5 


5 


39522 


166656 


3 


5 


4 


4332 


89280 


2 


5 


3 


14 


10416 


1 


3 


3 


252 


252 



code. This means that the [[7, 1, 5; 2]] and [[7, 1, 5; 3]] EAQEC codes cannot be obtained by 
the construction of [10]. Therefore, the [[7, 1, 5; 2]] and [[7, 1, 5; 3]] EAQEC codes are new. 

In addition, all the [[7, 1,5; 2]] EAQEC codes we found are degenerate codes, for some 
simplified stabilizer generators are of weight 4 from the check matrix. For example, a sim- 
plified check matrix and its simplified logical matrix of a [[7, 1, 5; 2]] EAQEC code are 



0000000 


1001011 


0000000 


1100101 


0000000 


0010111 


1001011 


0000000 


1100101 


0000000 


0010111 


0000000 


1000011 


0100011 


1101000 


0010010 



1001011 0100011 
1101000 1001011 



with T = {1,4}. On the other hand, all the [[7, 1,7; 6]] EAQEC codes are nondegener- 
ate codes, while the [[7, 1,5; 3]], [[7, 1, 5; 4]], [[7, 1, 5; 5]] EAQEC codes can be degenerate or 
nondegenerate. 



□ 
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Example 5. The check matrix of Shor's [[9, 1, 3]] code is 



000000000 
000000000 
000000000 
000000000 
000000000 
000000000 
111111000 
000111111 



110000000 
011000000 
000110000 
000011000 
000000110 
000000011 
000000000 
000000000 



As can be seen in Table HI the parameters [[9, 1, 9; 8]], [[9, 1, 7; 5]] and [[9, 1, 7; 4]] achieve the 
quantum singleton bound for EAQECC ( ITTj) and are not equivalent to any standard quantum 
stabilizer code. A competing EAQEC code with n = 9 and d = 7 by the construction of 
[TU] is a [[9, 1, 7; 6]] EAQEC code, obtained from a [9, 2, 7] linear quaternary code in Grassl's 
table. Therefore, the [[9, 1, 7; 5]] and [[9, 1, 7; 4]] EAQEC codes are new. All the [[9, 1, 5; 2]], 
[[9, 1,5; 3]], [[9, 1,7; 4]], [[9, 1, 7; 5]] and [[9, 1,7; 6]] codes are degenerate codes, and all the 
[[9, 1, 9; 8]] codes are nondegenerate codes, while the [[9, 1, 7; 7]] codes can be degenerate or 
nondegenerate. 



Table 4: Optimization over Shor's [[9, 1, 3] 



c 


dopt 


dstd 


N opt 


2' 2ck N(r, c) 


8 


9 


7 


256 


65536 


7 


7 


6 


330624 


4.17 x 10 6 


6 


7 


6 


278904 


4.42 x 10 7 


5 


7 


6 


17748 


9.94 x 10 7 


4 


7 


5 


132 


5.14 x 10 7 


3 


5 


5 


69777 


6.21 x 10 6 


2 


5 


5 


201 


1.72 x 10 5 



code 



□ 
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4 A Random Search Algorithm for Encoding Optimiza- 
tion 



It is easy to check that 

2 C (^ C ) < N(r,c) < Q2 c(r ~ c) . 

Hence 

A complete encoding optimization procedure for a [[n, k, d]} standard code becomes impossi- 
ble when n + k becomes large. Hence we consider a random search algorithm for the encoding 
optimization procedure in this section. 

We define the weight enumerator of the set N(S') — Sj as f(x) = Y17=o a i x% i where a; is 
the number of elements of weight i in N(S') — Sj. Suppose H' is a simplified check matrix 
of the EAQEC code corresponding to S'. Then we define a merit function m : H' i— >• M by 

b 

m(H') = y^a^, 

8=1 

where b = [ ra+c ~ fc ~ 2 j is an upper bound on the minimum distance for the given parameters 
n, k, c. This merit function was chosen to try to maximize the number of correctable errors in 
the set of likely errors for a typical memoryless channel; but of course, other merit functions 
could be chosen. The random search algorithm for encoding optimization is as follows: 

1) Apply Wilde's encoding circuit algorithm to a given check matrix H to obtain the en- 
coding operator Ue- 

2) Apply Ue to the simplified check matrix ( Tl3l) and the simplified logical matrix ( Tl6i) . and 
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10 

11 

12 



denote the new matrices by H' Q and L' Q , respectively. 
Compute the merit function m(H'^). 
Randomly choose a set T = {ii, • • • ,t c }. 

Randomly generate the matrices A and B according to Theorem |5] 

Apply a selection operator based on A and B, and obtain a new simplified logical matrix 
L'. 

Randomly generate a matrix My according to Theorem [3j 

Apply the unitary row operator defined by My to obtain a new simplified logical matrix 
L' . 

Compute the minimum distance d of the EAQEC code defined by H' and V . 
Compute the merit function m(H'). 
If m(H') < m(H' Q ), set H' Q < — H' . 

If d = b = [ n+c ~ fc ~ 2 j , or the search has repeated more than a maximum number of times, 
stop. 



13) Else, Go to step 4). 

Since the minimum distance might not be the best measure of the error-correcting ability of 
a quantum error-correcting code, a different merit function can be adopted. Here we simply 
chose the merit function m() to encourage EAQEC codes with fewer low-weight elements in 
N(S') — Sj. Some examples of this random search algorithm follow: 

Example 6. Applying the random search algorithm to the [[7,1,3]] quantum BCH code 
and Shor's [[9, 1, 3]] code, we obtain the same results as those in Table [3] and Table H] by 
the complete encoding optimization procedure. We list in Table E] and Table El the average 
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number of loops, N avg , to obtain these parameters (which we estimate by performing the 
algorithm 10 4 times for each c). Note that N avg is close to 2 2ck N(r,c)/N opt for the [[7, 1,3]] 
quantum BCH code, except for the case c = 2, where N avg is about half of 2 2ck N(r, c)/N opt . 
For Shor's [[9, 1,3]] code, N avg is much less than 2 2ck N(r,c)/N opt in most cases. 



Table 5: Average number of trials of the random search algorithm on the [[7, 1, 3]] quantum 
BCH code 



c 


dopt 


2 2ck N(r,c)/N opt 




6 


7 


113.78 


109.37 


5 


5 


2.02 


2.02 


4 


5 


4.22 


4.13 


3 


5 


20.61 


18.63 


2 


5 


744 


392.57 



Table 6: Average number of trials of the random search algorithm on Shor's [[9, 1, 3]] code 



c 


dopt 


2 2ck N(r,c)/N opt 


N a vg 


8 


9 


256 


259.97 


7 


7 


12.61 


12.47 


6 


7 


158.48 


119.25 


5 


7 


5600.63 


1188.08 


4 


7 


389393.94 


7577.95 


3 


5 


88.99 


18.48 


2 


5 


855.72 


52.46 



□ 



Example 7. We applied the random search algorithm to a standard [[15, 7, 3]] quantum BCH 
code and the results are shown in Table [71 Note that we did not find the fully optimized 



Table 7: Optimization over a [ 



c 


dopt 


d s td 


2 2ck N(r, c) 


8 


6 


5-6 


5.19 x 10 33 


7 


5 


5-6 


8.08 x 10 31 


6 


5 


5-6 


2.08 x 10 29 


5 


4 


4-5 


1.14 x 10 26 


4 


4 


4-5 


1.44 x 10 22 


7 


3 


4 


4.27 x 10 17 



15, 7, 3]] Quantum BCH code 



parameters in this case, since the complexity is very high. However, compared with the 

[[15, 3, 5; 4]] EAQEC code, obtained by the construction of ( fT9l) from a [15,7,5] classical 
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BCH code, the [[15, 7, 5; 7]] and the [[15, 7, 5; 6]] EAQEC codes have 4 more information 
qubits at the cost of 3 and 2 more ebits, respectively. The [[15, 7, 6; 8]] EAQEC code has 
4 more information qubits and a higher minimum distance at the cost of 4 more ebits. In 
addition, the [[15, 7, 6; 8]] EAQEC code is not equivalent to any known standard quantum 
stabilizer code. 

On the other hand, the classical linear quaternary [15, 9, 5] code and [15, 8, 6] code 
in Grassl's table can be used to construct a [[15, 9, 5; 6]] EAQEC code and a [[15, 9, 6; 8]] 
EAQEC code by the construction of [TU], respectively. These two codes are better than 
the [[15, 7, 6; 8]] EAQEC code we obtained. This may be because these codes were not fully 
optimized, but BCH codes in any case need not give the best possible EAQEC codes, even 
using the encoding optimization procedure. 

□ 

Example 8. We applied the random search algorithm to Gottesman's [[8, 3, 3]] code jl] and 
the [[13, 1, 5]] quantum QR code [21 HE], and the results are shown in Table and Table [9J 
respectively. By the construction of [TU], the [8,3,5], [13,3,9], [13,4,8], [13,5,7] classical 
linear quaternary codes in Grassl's Table can be transformed to [[8, 2, 5; 4]], [[13, 3, 9; 10]], 
[[13, 0,8; 5], [[13, 1,7; 4]] EAQEC codes, respectively. Hence the [[8, 3, 5; 5]], [[13, 1, 11; 11]], 
[[13,1,11; 10]], [[13, 1,9; 9]], [[13, 1,9; 8]] EAQEC codes are new, and are not equivalent to 
any standard quantum stabilizer code. □ 



Table 8: Optimization over Gottesman's 



c 


dopt 


dstd 


2 2cA W(r, c) 


5 


5 


4 


1.07 x 10 9 


4 


4 


4 


5.20 x 10 7 


3 


4 


3 


4.06 x 10 7 


2 


3 


3 


6.34 x 10 6 



[8, 3, 3]] code 
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Table 9: Optimization over the [[13, 1, 5]] quantum QR code 



c 


dopt 


d s td 


2**jV(r, c) 


12 


13 


9 


1.68 x 10 7 


11 


11 


8-9 


1.71 x 10 10 


10 


11 


7-9 


2.92 x 10 12 


9 


9 


7-8 


1.07 x 10 14 


8 


9 


7 


9.12 x 10 14 


7 


7 


7 


1.87 x 10 15 


6 


7 


7 


9.44 x 10 14 


5 


7 


7 


1.17 x 10 14 


4 


7 


7 


3.56 x 10 12 



5 Circulant Construction of EAQEC codes 



The construction in Theorem [2] and the construction (fP9"j) resemble the CSS construction 
[7J [8] , while standard quantum stabilizer codes with a high minimum distance are usually 
non-CSS codes. In this section we give a different construction of EAQEC codes, with a 
simplified check matrix that differs from the simplified check matrix in construction ( Ti~9l . 
We construct the simplified check matrix directly, rather than starting from a classical binary 
code. 

Let H' be a r x 2n simplified check matrix cyclicly generated by a binary 2n— tuple 

a = ao«l • • • 0,2n-2Q<2n-l'- 



On • • • «2n-l 

On+1 «n+2 • • • «0 



^r+n— 1 O'r+n ■ ■ ■ &r—2 

If the rank of H' is exactly r, then c = |rank(i7'Ai?') and if' defines an [[n, n + c — r, <i; c]] 
EAQEC code for some minimum distance d. For example, a [[6, 1,4; 1]] code is constructed 



H'=[H' X \H' Z ] 



do a\ . . . a n -i 



ai a 2 



a,, — i a r ... a r + n _2 
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by a = 001110101110 with the simplified check matrix 



001110 


101110 


000111 


010111 


100011 


101011 


110001 


110101 


111000 


111010 


011100 


011101 



We call this the circulant construction of EAQEC codes, which is used for standard stabilizer 
codes in [15] . 

We examined the simplified check matrices cyclicly generated by every possible binary 
2n— tuple a by computer for n = 4, ■ • • , 10 and r < 2(n — 1). Parameters of EAQEC 
codes not equivalent to any standard quantum stabilizer codes are listed in Table [TUJ The 
parameters [[4, 0,4; 2]], [[4,1,3;1]], [[5, 0,4; 2]], [[5, 1,4; 3]], [[5, 1,5; 4]], [[6, 0,6; 4]], [[6,2,3; 1]], 
[[6, 2, 4; 3]], [[6,0,4; 1]], [[6, 1,5; 4]], [[7, 4, 3; 2]], [[7, 1,6; 5]], [[7, 1,7; 6]], [[8, 5, 3; 2]], [[9, 1,9; 8]], 
and [[10, 0, 10; 8]] also saturate the quantum singleton bound ffTTl) . 

Table 10: Parameters of [[n, k, d; c]] EAQEC codes not equivalent to any standard [[n + c, k}} 
codes. 



n 


[[n,k,d;c]] 


4 


[[4, 0,4; 2]], [[4,1,3; 1]] 


5 


[[5, 1,5; 4]], [[5, 1,4; 3]], [[5, 1,4; 2]], [[5, 0,4; 2]], [[5, 2, 3; 2]] 


6 


[[6, 0,6; 4]], [[6, 1,5; 4]], [[6, 1,4; 3]], [[6, 2, 4; 3]], [[6,0,4; 1]], [[6,2,3;1]] 


7 


[[7, 1,7; 6]], [[7, 2, 5; 5]], [[7, 0,6; 4]], [[7, 3, 4; 4]], [[7, 1, 4; 2]], [[7, 3, 4; 3]], [[7, 4, 3; 2]] 


8 


[[8, 0,8; 6]], [[8, 1,6; 6]], [[8, 0,6; 5]], [[8, 2, 6; 6]], [[8, 1,6; 5]], [[8, 0,6; 4]], 
[[8, 3, 5; 5]], [[8, 2, 5; 4]], [[8,1,4; 1]], [[8, 3, 4; 3]], [[8, 5, 3; 2]], 


9 


[[9, 1,9; 8]], [[9, 0,7; 6]], [[9, 1,7; 6]], [[9, 1,7; 7]], [[9, 2, 6; 6]], [[9, 1,6; 5]], 
[[9, 0,6; 4]], [[9, 1,6; 6]], [[9, 2, 5; 4]], [[9,5,3;1]], 


10 


[[10, 0,10; 8]], [[10, 1,8; 8]], [[10, 0, 8; 7]], [[10, 0, 8; 6]], [[10, 0, 7; 5]], [[10, 1, 7; 6]], 
[[10, 2, 7; 7]], [[10, 1,6; 5]], [[10, 3, 6; 7]], [[10, 0, 6; 3]], [[10, 3, 6; 6]], [[10, 2, 6; 5]], 
[[10, 1,6; 4]], [[10, 4, 5; 5]], [[10, 2, 5; 2]], [[10, 4, 5; 4]], [[10, 2, 5; 3]], 
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6 Discussion 



This paper has studied how entanglement can be used to increase the minimum distance 
of quantum error-correcting codes. We demonstrated the encoding optimization procedure 
for EAQEC codes obtained by adding ebits to standard quantum stabilizer codes. The 
four types of unitary row operators play an important role in this encoding optimization 
procedure, and also help to clarify the properties of EAQEC codes and their relationship to 
standard codes. Some applications of the encoding optimization procedure were found to 
have promising results: we found [[7, 1,5; 2]] and [[7, 1,5; 3]] EAQEC codes from quantum 
BCH codes, and [[9, 1, 7; 4]] and [[9, 1, 7; 5]] EAQEC codes from Shor's 9-qubit code, together 
with a family of [[n, 1, n; n — 1]] EAQEC codes for n odd, all of which achieve the quantum 
singleton bound. Several EAQEC codes found by this encoding optimization procedure are 
also degenerate codes. This procedure serves as an EAQEC code construction method for 
given parameters n, k, c. 

The encoding optimization procedure has very high complexity. However, it might 
be useful to further investigate it for specific families of codes that have special algebraic 
structures, such as quantum BCH codes and quantum Reed-Muller codes. For example, if 
we added ebits to a quantum Reed-Muller code (by the CSS construction with a classical 
RM(r, m) code [2]), we found that its simplified logical matrix V is, like a check matrix of 
a CSS code, of the form 

O L' z 
L'x O 

By examining the symplectic relations, we found that the classical code, generated by L' z 
or L' x , together with the RM(r,m) code, generate a subcode of RM(r',m), which con- 
tains more lower-weight codewords. The minimum distance of the EAQEC code might be 
increased approximately by a factor | similar to the method used in [19] or [18]. From a 
[[16, 6, 4]] quantum Reed-Muller code, we obtained a [[16, 6, 6; 10]] EAQEC code by construct- 
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ing the matrix My such that the above argument could be applied. However, we obtained 
a [[16, 6, 7; 10]] EAQEC code by applying the random search algorithm for the encoding op- 
timization procedure. How to directly construct a matrix My that leads to EAQEC codes 
with high minimum distances is a subject of ongoing research. 

When the complexity becomes large, it is almost impossible to optimize over all 2 2ck N(r, c) 
encoding operators. The random search algorithm seems to be the only method to achieve 
good (but suboptimal) results for EAQEC codes. For different parameters n, k, c, the merit 
function should be carefully chosen. The best choice of merit function for a given applica- 
tion is also a subject of future work. A search algorithm for specific EAQEC codes could be 
developed. While the encoding optimization procedure in this paper applies to a standard 
quantum stabilizer code, it is possible to construct a similar encoding optimization algorithm 
for adding ebits to other EAQEC codes that have ancilla qubits which are not ebits. Much 
work remains to be done in finding the best possible EAQEC codes for different applications. 
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